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Aims. We quantify the contributions of 24 (j,m galaxies to the Far-Infrared (FIR) Background at 70 and 160 /im. We provide 
new estimates of the Cosmic Infrared Background (CIB), and compare it with the Cosmic Optical Background (COB). 
Methods. Using Spitzer data at 24, 70 and 160 fim in three deep fields, we stacked more than 19000 MIPS 24 /im sources with 

■ S24 > 60/tJy at 70 and 160 /xm, and measured the resulting FIR flux densities. 

"**••*»• Results. This method allows a gain up to one order of magnitude in depth in the FIR. We find that the Mid-Infrared (MIR) 
24 ^m selected sources contribute to more than 70% of the Cosmic Infrared Background (CIB) at 70 and 160 ^m. This is the 
first direct measurement of the contribution of MIR-selected galaxies to the FIR CIB. Galaxies contributing the most to the 
Q ' total CIB are thus z ~ 1 luminous infrared galaxies, which have intermediate stellar masses. We estimate that the CIB will be 
?—j \ resolved at 0.9 mjy at 70 and 3 mjy at 160 /im. By combining the extrapolation of the 24 fim source counts below analysis, we 
CZJ . obtain lower limits of 7.1 ± 1.0 and 13.4 ± 1.7 nW m~ 2 sr" 1 for the CIB at 70 and 160 /im, respectively. 

Conclusions. The MIPS surveys have resolved more than three quarters of the MIR and FIR CIB. By carefully integrating 
the Extragalactic Background Light (EBL) SED, we also find that the CIB has the same brightness as the COB, around 

■ 24 nW m~ 2 sr -1 . The EBL is produced on average by 115 infrared photons for one visible photon. Finally, the galaxy formation 
and evolution processes emitted a brightness equivalent to 5% of the primordial electromagnetic background (CMB). 

Key words. Cosmology: observations - Cosmology: Diffuse Radiation - Galaxies: Evolution, Starburst, Infrared 



1. Introduction 

The Cosmic Infrared Background (CIB) is the relic emis- 
sion at wavelengths larger than a few microns of the for- 
mation and evolution of the galaxies of all types, including 
Active Galac t ic Nu c lei (AGN) and star- forming; systems 
llPuget et all Il996l iHauser et all Il998l iLaeache et all 



19991 ICispert et all 1200(1 lHauser fc DwekL l200lt 
KashlinskvL 120051) . Characterizing the statistical behavior 
of galaxies responsible for the CIB - such as the number 
counts, redshift distribution, mean Spectral Energy 
Distribution (SED), luminosity function, clustering - 
and their physical properties - such as the roles of 
star-forming vs accreting systems, the density of star 
formation, and the number of very hot stars - has thu s 
been an important goal l|Partridge fc Peeblesl Il967jl . 
The SED of the CIB peaks near 150 fim. It accounts for 
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roughly half of the total energy in the optical /infrared 
Extra galactic Background Light (EBL) (Hauscr & Dwck, 
2QQ1) , although still with some uncertainty IjWrighti 



2004: lAliaronian et all Eo05h . Since locally the infrared 



output of galaxies is only a third of the optical one 
llSoifer fc Neugebauerl Il99ll) . there must have been a 
strong evolution of galaxy properties towards enhanced 
Far-Infrared (FIR) output in the past. Understanding 
this evolution requires interpretation of cosmological 
surveys conducted not only in the infrared and submil- 
limeter spectral rang es, but also at other wavelengths 
l|Lagache et all 12005^ . 

The cryogenic infrared space missions IRAS (Infrared 
Astronomical Satellite) and ISO (Infrared Space 
Observatory) provided us with valuable insights to 
the IR-domi nated galaxies in the Mid-Infrared (MIR) 
and FIR llSanders fc Mirabel I"l996t iGenzel fc Cesarskvl 
l2000t iDolel 120031 lElbazl 12003 lLagache et all l2005l for 
reviews). ISO MIR surveys were able to resolve a signifi- 
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Fig. 1. Sensitivity to the bolometric luminosity (and star 
formation rate, assuming star forming galaxies) of var- 
ious infrared and submillimeter experiments. Detections 
of at least 10 sources in the surveys can be made in 
the areas above the curves. We assumed the scenario of 
a typical deep survey. ISOCAM 15 fim {S v > 250/iJy, 
2 Sq. Deg.); ISOPHOT 170 fim (S„ > 180 mJy, 5 
Sq. Deg.); Spitzer/MIPS 24 fim (S v > 80/iJy, 5 Sq. 
Deg.); Spitzer/MIPS 70 fim (S v > 25 mJy, 5 Sq. Deg.); 
Spitzer/MIPS 160 fim (S v > 50 mJy, 5 Sq. Deg.); SCUBA 
850 fim (S v > 1 mJy, 1 Sq. Deg.); This plot makes use of 
the Lagache et al. (2004) model and their starburst SED 
for the conversion to Lboi- At z ~ 1, MIPS detects only 
ULIRGs in the FIR, and detects LIRGs in the MIR. The 
stacking analysis allows to gain an order of magnitude and 
to probe LIRGs in the FIR. 



cant fraction of the 15 /im CIB ijElbaz et allll999l counts 
close to conver g ence). Using model SEPs of galaxies 
llCharv fc Elbad. l200lt IXu et all l200lt lLagache et all 
120031 for instance), the contribution of MIR-selected 
galaxies to the peak of the CIB (a round 140 to 170 /im) 
can be inferred. lElbaz et alJ l|2002|) derived that 64 ± 38% 
(16 ± 5 over 25 ± 7 nW m" 2 sr" 1 ) of the 140 fim 
background is due to ISOCAM 15 fim galaxies, whose 

median redshift is z ~ 0.8. 

The Spitzer Observatory I Werner et all l2004h is per- 
forming much deeper and wider-area surveys, in par- 
ticular at 24, 70 to 160 fim using the Multiband 
Imag ing Photometer for Spitzer (MIPS) ijRieke et all 
|2004[) . However, because of the limited angular resolution 
( "smoothing" the high spatial frequency signal in the FIR 
map s), deep MIPS 70 and 160 fim maps are confusion lim- 
ited ijDole et alll2003ll2004bh - the source surface density 
corresponds to 20 beams per source or less e.g. in the GTO 
fields. The FIR images do not allow us to directly probe 
the same galaxy population as that detected at 24 fim, 
where the extragalactic source confusion is less important. 
Figure shows the typical sensitivity of MIPS surveys to 
the bolometric luminosity of galaxies as a f unction of red- 
shift, using the modeled starburst SED of lLagache et al.1 



( 2004). At a redshift z ~ 1, MIPS FIR surveys are sensi- 
tive to ultraluminous IR galaxies (ULIRGs, L > 10 12 L Q ) 
where MIPS 24 fim surveys can probe luminous IR galax- 
ies (LIRGs, L > 1O 11 L0). It is therefor impossible to de- 
rive MIR and FIR SEDs of individual LIRGs at z ~ 1 and 
above. 

MIPS can detect high redshift sources at 24 fim: 
about 25 to 30% of the population of galaxies lie 
at z > 1.5 . at faint flux densities (down to fe w 
tens of uJv) IjLe Floc'h et all Eooi lEgami et all l20(L 
Lonsdale et all l2004 ICIwv et all l2004 iHonck et a,l 



2005iTperez-Gonzalez et all l2005t ICanuti et alT E 
Papovich et alJ 1)20041) showed that MIPS surveys resolve 
about 70% of the 24 fim IR galaxy CI B for S9.4 > 60 mJv . 
In comparison, MIPS 70 and 160 fim l|Pole et all 12004^ 
surveys detect only about 20% and less than 10% of the 
CIB at 70 and 160 fim, respectively. Programs of very 
deep 70 fim imaging on small fields (e.g. D. Frayer, private 
communication) are likely to resolve a larger fraction, but 
due to confusion noise the FIR CIB will still not be sig- 
nificantly resolved into individual sources, while the MIR 
CIB at 24 fim is well resolved. 

In this paper, we use a stacking analysis method that 
takes advantage of the good sensitivity of the MIPS 24 fim 
MIR channel, to fill the sensitivity gap between the MIR 
and the FIR surveys. By stacking the FIR data at the 
locations of MIR sources, we statistically investigate the 
FIR properties of 24 /im-selected galaxies. In particular, 
we quantify the contribution of the 24 fim resolved galax- 
ies to the 70 and 160 fim background, put strong lower 
limits to the CIB, and give new estimates of the 70 and 
160 fim background. 

Throughout this paper, we adopt a cosmology with 
h = 0.65, Q.m = 0.3 and £!a = 0.7. The surface bright- 
nesses (e.g. of the CIB) are usually expressed in units of 
MJy/sr or nW m -2 sr -1 . For a given frequency v in GHz 
and wavelength A in microns, the conversion between the 
two is given by: 



InW m ~ sr 



100 



v/GHz 



MJy/sr ■ 



X/ fim 
3000 



MJy/sr (1) 



2. Data and Sample 



The data for our analysis are from the Spitzer MIPS 
Guaranteed Time Observations (GTO) cosmological sur- 
veys performed in three fields: the Chandra Deep Field 
South (CDFS), the Hubble Deep Field North (HDFN) 
and the Lockman H ole (LH). The MIPS o bservations at 
24 fim are detailed in Paj | ) ovich et al 1 1|2004 and at 70 and 
160 fim in lPole et"aT. 1 2004^ . Each field covers about 0.4 
square degrees, and the integration times per sky pixel 
are 1200s, 600s and 120s at 24, 70 and 160 fim, respec- 
tively. The data wer e reduced and mosa icked using the 
Pata Analysis Tool ((Gordon et all l2005h . We make use 
of a recent new analysis of the calibration by the instru- 
ment team and the instrument support team that will soon 
be adopted officially. The uncertainty is now 4%, 7% and 
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12% at 24, 70 and 160 /im, respectively, and the calibra- 
tion level has been changed by less than 10% compared 
with the previous determinations (See the Spitzer Science 
Center calibration pages 1 ). 

IPanovich et all l|2004h showed that the 80% complete- 
ness level at 24 /im in the GTO deep fields is rea ched 
at 5*24 = 80/ i Jy. Ne vertheless, IPapovich et alJ l|2004j) and 
ICharv et al.l l|2004|) show that the detection of very faint 
24 (im sources, down to 5*24 ~ 30/iJy, is possible, but with 
a increased photometric uncertainties and reduced com- 
pleteness (to lower than 5% at the GTO depth and 20% 
at the GOODS depth) . 

bole et alJ (|2004al) showed that at 70 and 160 /im 
sources can be safely ex tracted down to 15 mJy and 
50 mJy, respectively. The iFraver et alJ l|2006() results go 
deeper. However, confusion limits the extraction of sources 
fainter than typically 5 6 /iJy at 24 /im, 3.2 mJy at 70 /im, 
and 36 mJy at 160 /im (Do le et al A priori infor- 

mation on the existence of a source deduced from shorter 
wavelength and less confusion-limited observations can ex- 
tend the reliable detection threshold below this nominal 
confusion limit. 

To implement this approach, we build a sample as fol- 
lows: 

• We select the central part of each field where all 3 MIPS 
wavelengths have a common sky coverage and maximum 
redundancy. This area covers 0.29 Sq. Deg in the CDFS 
and LH, and 0.27 Sq. Deg. in HDFN, for a total of 0.85 
Sq. Deg. in these three fields. 

• In these selected areas, we identify every MIPS 24 /im 
source with S24 > 60/iJy. This flux density limit cor - 
responds to 50% completeness (jPapovich et all |2004). 
There are 6543 galaxies above this limit in CDFS, 6039 
in the HDFN, and 6599 in the LH. The total number of 
sources considered in the three fields is thus 19181. 

To analyze this sample, we proceed as follows: 

• In each field, we sort the 24 /im sources by decreasing 
flux density 6*24. 

• We put the sources in 20 bins of flux density for 
S*24 > 60/iJy. These bins have equal logarithmic width 
AS24/S24 ~ 0.15, except for the bin corresponding to 
the brightest flux, which includes all the bright sources 
(0.92mJy to Uy). 

• We correct the average flux obtained by stacking each 
i5*24 bin for incomple teness, following the correction of 
IPapovich et all l|2004|) (their figure 1). Since the bins be- 
tween 60 and 80 /iJy are complete to the 50-80% level, 
only the weakest fluxes bins are significantly corrected. 

3. Stacking Analysis 

The process of stacking the sources based on the 24 /im 
detections allows us to measure more of the total contri- 
bution of 70 and 160 /im sources to the CIB. 

1 http:/ /ssc. spitzer. caltech.edu/mips/calib/conversion. html 




Fig. 2. Images at 24, 70 and 160 /im (left to right) of 
stacked sources in the brightest bin of 24 /im flux den- 
sity, with a random position offset added before each sum. 
Color coding: dark is high flux, light is low flux. This al- 
lows us to check that the stacking method does not in- 
troduce any artifacts, i.e. a false source detection in the 
center. 

3.1. Processing 

At 24 /im, the detector pixel size is 2.5 arcseconds, the 
FWHM of the point spread function (PSF) is 6 arcsec- 
onds, and the plate scale of the mosaic is chosen to be 1.25 
arcseconds. At 70 /im, the detector pixel size is 9.9 arc- 
seconds, the FWHM of the PSF is 18 arcseconds, and the 
plate scale of the mosaic is chosen to be 4.5 arcseconds. 
At 160 /im the detector pixel size is 18 arcseconds, the 
FWHM of the PSF is 40 arcseconds, and th e plate scale 
of th e mosaic is chosen to be 18 arcseconds l|Rieke et all 
l2004t iGordon et all Eooi for more details). The 70 and 
160 /im mosaics have been resampled to the scale of the 
24 /im mosaic (1.25 arcseconds per pixel) using a bilinear 
interpolation. This last step greatly facilitates the weight 
management of the three maps, since each has different 
coverage, and it allows easy extraction of the signal at the 
three wavelengths for the same sky position. 

For each S24 flux density bin we select every 24 /im 
source, extract a square image about 440 arcseconds on 
a side centered on the source, and store it. We proceed 
similarly on the mosaics at 70 and 160 /im, extracting 
images at the position of each 24 /im source regardless 
of the presence of a detected FIR source. The products 
at this stage are thus three cubes of data at 24, 70 and 
160 /im with the same dimensions (same number of source 
images and same box size) for each of the 24 /im flux 
density bins. 

We then add the images in each cube at each wave- 
length, to generate a stacked image of sources at 24, 70 
and 160 /im for a given S 2 4 flux bin. This operation is a 
simple sum, without any outlier rejection. When stacking, 
we rotate each image by +7r/2 with respect to the previous 
one (and so on), to cancel out the large-scale background 
gradients such as the prominent zodiacal background at 
24 /im. This processing is done both in each field sepa- 
rately as well as using all the data at once. Unless oth- 
erwise stated, we use the stacked data of all the fields 
together in the rest of this paper. We checked that no 
significant signal was detected when we added a random 
or systematic artificial offset to each 24 /im position and 
then performed exactly the same sub-image extraction 
and stacking as we did for the real 24 /im source list. 
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Fig. 3. Images at 24 //m of stacked sources in bins of 
24 /im flux density (SW > 60^Jy) in the three MIPS 
GTO Fields: CDFS, HDFN, and LH covering about 0.85 
Sq. Deg. A total number of 19181 sources has been used. 
The number of sources used in the sum in each S24 bin is 
reported. Each image has 350 x 350 pixels of 1.25 arcsec, 
thus covering about 7.3 x 7.3 Sq. Arcmin. Since no outlier 
rejection has been made, other sources can be seen in the 
surroundings of the stacked sources 



Figure |21 shows the results of stacking the sources (only 
from the brightest bin), with a random position offset 
added prior to the sum. No source appears in the center, 
as expected. This guarantees that the stacking method 
does not introduce an artifact that mimics a source. 

Since the stacking analysis aims at statistically detect- 
ing faint unresolved sources at 70 and 160 /im, in princi- 
ple there is no need to also stack data at 24 /im, where 
all sources are resolved. However, doing so allows us to 
double-check the method, since we know by design what 
the stacked photometry should be. 

3.2. Stacked Images and Photometry 

The final stacked images at 24, 70 and 160 /im as a func- 
tion of the 24 /im flux density S24 are presented in fig- 
ures |31 0J and respectively. We report also in these fig- 
ures the number of sources stacked in each of the S24 bins. 
The figures show clear detections of stacked sources at 70 
and 160 /im for every S24 bin, even the faintest corre- 
sponding to 60 < 5 2 4 < 69 /iJy. Given the surface den- 
sity o f the 24 urn sou r ces at 60 fiJy of (9.6 ± 0.04) x 10 7 
sr -1 fPapovic h et al ., 2004). this translates to 1.04 and 





Fig. 4. Images at 70 /im of stacked sources in bins of 
24 /zm flux density. See caption of Figure 01 for details. 



0.2 beams per so urce at respectively 70 and 160 /im 
l|Dole et al.L l2003. using beams from their Table 1). This 
is well beyond the confusion limits at these FIR wave- 
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Fig. 5. Images at 160 /im of stacked sources in bins of 
24 /im flux density. See caption of Figure [3] for details. 
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Fig. 6. Normalized radial profiles of the stacked images. 
The crosses represent the data, and the solid line the em- 
pirical PSF. The vertical dotted lines show the radii of 
the aperture used for photometry. From top to bottom: 
faintest S24 bin at 70 /xm; Brightest 6*24 bin at 70 /xm; 
Faintest S24 bin at 160 /xm; Brightest S24 bin at 160 /xm . 

lengths l|Dole et all l2004bl) . This statistical detection of 
FIR sources already demonstrates the great potential of 
this technique to probe FIR galaxies down to levels be- 
low the confusion, thanks to the excellent quality of the 
pointing and the stability of the effective PSF (see below). 

We check that the radial profile of the stacked sources 
is in agreement with the PSF profile, at each wavelength 
and for each flux bin. We show in Figure two profiles 
at each wavelength corresponding to the extreme cases: 
the brightest and faintest S24 flux bins. We used both 
the empirical PSF (from bright sources) and the mod- 
eled STinyTim MIP S PSF (jlOisflll 993llRieke et allEooH 
iGordon et all 120051) . At large S24, the stacked radial pro- 
files at 70 and 160 (im (bottom plots in Figure^ agree well 
with the PSF in the central part. At the faintest fluxes (top 
plots), the agreement is good down to about 10% of the 
peak brightness. Since the stacked images visually repre- 
sent the 2-dimensional correlation function of galaxies, the 
potential presence of many neighboring sources at small 




0.1 1.0 
Flux of 24 um Sources [mjy] 

Fig. 7. Top: Confidence Level of the detections at 70 
(dash) and 160 /xm (solid) on stacked images (Fig. ^ and 
EJ, as a function of the S24 bin. Middle: Signal-to-Noise 
ratio, as computed from a Gaussian fit to the flux distri- 
bution measured on about 2000 positions; this S/N is not 
relevant at low flux (see Sect. I3.2j> . Bottom: average flux 
in mjy per stacked source. Note that a different number 
of sources have been stacked in each bin. 



scales (source clustering) might have widened the radial 
profile, which is not observed; thus source clustering does 
not contribute significantly to the noise budget. 

We measure the flux density of the stacked sources 
with aperture photometry and correct for aperture size. 
The radii of the apertures and reference annulus are, in 
arcseconds: (r a p er , ri n t, r ex t ) = (12.2, 17, 24), (30, 49, 
79) and (54, 90, 126) for 24, 70 and 160 /xm, respectively. 
These radii correspond to approximatively 3, 5 and 7 times 
the FWHM in the FIR and 2, 3, 4 times the FWHM in 
the MIR. We measured the noise in each image by us- 
ing about 2000 measurements on random positions. We 
compute the confidence level (C.L.) of each detection (top 
of Figure 01 using the cumulative distribution of the noise 
measurements. The deviation from 100% of the C.L. is the 
probability that the noise creates a spurious source. For 
the faintest bin, the C.L. is around 80%, and it rises to 97% 
for the next four bins, and stays at 100% for the brighter 
S24 bins. We fitted a Gaussian function to the distribution 
of noise to get the standard deviation in order to estimate 
the S/N ratio. This method works for the brighter bins 
(middle panel in Figure EJ), where the flux distribution 
is indeed nearly a Gaussian distribution. In this range, 
the S/N values have a median of 8 at 70 /xm and 7 at 
160 /im. In the three faintest bins, the noise distribution is 
not Gaussian, because of the presence of slightly brighter 
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Fig. 8. Contributions to the CIB: brightness of stacked 
sources at 160 and 70 fim per logarithmic bin dl ^g 2 ^ 
(solid symbols), as a function of the 24 /zm flux, in all 
three MIPS GTO Fields. A completeness correction has 
been applied. The highest flux bin goes up to Uy. Open 
symbols: published differential source counts multiplied by 
S^' 5 and a color ratio of bright galaxies (160/24=60 and 
70/24= 20, cf the br i ghtest b in in the lower panel of Fig.pT ); 
Square: iDole et all l(2004ah : Circle: iFraver et all (Hl)06). 
There is a good agreement between the source counts, the 
brightest stacked bin, and the fainter stacked bins. 



sources; the Gaussian fit is therefore not relevant and we 
opt for the C.L. technique. The bottom plot in Figure 7 
shows the average FIR flux per stacked galaxy. A set of 
~100/iJy MIR-selected galaxies would have a typical av- 
erage FIR flux of ~ 0.5 and 3 mjy if taken individually at 
70 and 160 /im, respectively. Since the confusi on limits are 
at abo ut 3 and 40 mjy at these wavelengths l)Dole et all 
2004b), the gain of the stacking analysis technique is one 
order of magnitude in flux compared to individual detec- 
tion. Finally, it is not necessary to remove the brightest 
sources for the goals of this paper, because we stack typi- 
cally 1000 to 2000 galaxies per flux bin, so their influence 
is negligible except maybe in the 3 faintest bins. 

The brightness of the stacked sources at 70 and 160 



per logarithmic flux density bin, or 



dB x 



as a func- 



dlog(S) ' 

tion of the 24 /im flux bin, is presented in figure |H1 B\ 
in MJy/sr is defined as the total stacked flux density di- 
vided by the survey area. Using a logarithmic flux density 
bin allows direct comparison of the contribution in energy 
of each bin to the CIB, and is directly related to the dif- 
ferential source counts with a scaling factor S~ 5 . In the 
range 100 < S24 < 300 //Jy, both contributions to the CIB 
present a maximum, which shows that the contributions 



have reached convergence. Converting 5Vo and 5x60 into 
5 2 4 using the color ratios of 9 and 30 (see Table [21 be- 
low), this means the FIR CIB will be mainly resolved at 
S70 - 0.1 x 9 = 0.9 mjy and at S 160 ~ 0.1 x 3 = 3 mjy 
We ha ve al so plotted the sou rce counts of IDole et all 
( 2004^) and IFraver et alJ l(2006h : we used the conversion 
to S24 as given by color ratios relevant for bright galaxies 
of 20 and 60 at 70 and 160 /im, measured on the very 
bright end of the bottom plot in Figure [7| Despite this 
simplifying assumption of a single color ratio, there is ex- 
cellent agreement between the brightness derived from the 
stacking analysis and the source counts. This plot can be 
used to constrain models of galaxy evolution. 

Sample variance plays a role in these results. To probe 
its effects, we split each of our three fields (CDFS, HDFN, 
LH) into four subficlds of about 250 square arcmin each, 
and performed an independent analysis on each of these 
twelve subfields. We obtain contributions varying in some 
cases by as much as a factor of two (peak-to-peak). For in- 
stance, computing the standard deviation of the distribu- 
tion of the cumulative 160 flux (the faintest points in 
Figure IT7T|l measured over these 12 subfields gives a = 0.3 
MJy/sr and a mean and median both of 0.53 MJy/sr. 
Rcnormalizing by the twelve sub-fields gives a = 0.09: the 
uncertainty induced by the Large Scale Structure varia- 
tions across the fields is of order 15%. 

From here on in this paper, our error budget takes 
into account: 1) the calibration uncertainties; 2) the pho- 
tometric uncertainty; 3) the large-scale structure (sample 
variance). 



4. Contributions of Mid-Infrared Galaxies to the 
Cosmic Infrared Background 

4.1. Value of the Cosmic Infrared Background 
Brightness 

To compute the fraction of background resolved with 
the stacking analysis of the MIPS data, we first need to 
review the measurements of the total CIB, in particular 
at 24, 70 and 160 /im. It should be remembered that the 
total cosmic background contains the contribution of all 
extragalactic sources but also more diffuse emissions, e., 
from dust in galaxy clusters l|Montier &: Giardl I200E 
Furthermore the extragalactic sources are expected to 
be mostly galaxies but it cannot be excluded that other 
lower luminosity sources, population III stars for instance, 
contribute significantly but will not be detected directly 
in the present deep surveys. 

Measuring the CIB directly by photometry is partic- 
ularly difficult because one needs 1) an absolute pho- 
tometer and 2) a proper estimate of the foreground. The 
two FIR channels of MIPS are not absolute photome- 
ters for the very extended spatial scales, since no inter- 
nal calibrated reference can be observed to calibrate abso- 
lutely the slow response. A better knowledge of the instru- 
ment in the future may allow a proper absolute calibra- 
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Fig. 9. Current measurements of the Extragalactic Background Light Spectral Energy Distribution from 0.1 /mi 
to 1 mm, showing the Cosmic Optical Background (COB, with A < 8 /im) and the Cosmic Infrared Background 
(CIB, with A > 8 /im). Black arrows re present lower limits. P urple arrows and lines represen t upper limits. The 
EB L observational constr aints come from: lEdelstein et alJj2000t) at 0.1 /im using V oyager UVS; iBrown et alJ l|200Cfl 
andlGardner et alJ ll2000h with HST/STI S [lower limi tsl: iMadau fc PozzettH (l2000f) with HST (incl NICMOS) and 
iThomoson I (^2003^ 1: Bernstein et all fe002l) corrected by iMattilal l|2003|) [filled circles] : iMatsumoto et alJ (|2005h between 
2.2 and 4 ;xm using the IRT S [thin plus]: lGoriian et all (|200C|) at 2.2 and 3.3 /im using DIRBE and Lzcfc: IWrighl] d200lJ) 



_4i_ 

andlCambresv et all l)200l|) at 1.25 and 2.2 /im using DIRBE and 2MASS [five branch star]; DIRBE values from lWrightl 
112001" from 1.25 to 240 /im [gray circles]; Spitzer IRAC 3.6, 4.5, 5.8 and 8.0 u m lower limits fr om number counts 



bv lFazio et~aH (|2004h : fluctuation analysis with IRAC bv lSavage fc Oliver! l|2005h (open trianglesi : ISchroedterl ll2005h 
using Very High Energy Blazars, 98% confidence uppe r limit [gray region]; H .E.S.S upper limit from lAharonian et alJ 
l(2005j) using P0.55 [solid line b etween 0.8 and 4 /i m] ; iRenault et alJ l|200ll) upper limits from 5 to 15 /im using the 
CAT in the 7-rays on Mkn501: lElbaz et all I l9991 lower limit at 15 /im using galaxy counts with ISOCAM; upper 
limit at 20 /im bvlStecker fc De^agerTifl997|l on Mkn421: lower limit from galaxy counts at 24 /tm with MIPS by 
IPapoyich et alJ l)2004|k an indirect evaluation at 60 /im using fluctuations in IRAS data from lMiville-Deschenes et all 
l|2002|) [open gray square]; lower limits at 70 and 160 p,m using galaxy coun ts with MIPS bv lDole et al.1 ll2004al): an 
estimate of the CIB at 100 /mi using CAT and g/flffff lRenault et all l|200j) [four branch starl : lLagache et al l (|200oT) 
at 100, 140 and 240 /im using DIRBE an d WHAM, updated in the present work [diamond] ; lHauser et all l)l998|) at 140 
and 240 ^tm using D IRBE [open square] : ISmail et alJ l)2002j) lower limit at 850 /tm using galaxy counts with SCUBA; 
lLagache et alJ l)2000|) spectrum between 200 /im and 1.2mm using FIRAS [solid line above 200 /im]. The IDL script 
to generate this figure is available on the web: http://www.ias.u-psud.fr/irgalaxies 



tion, using its "Total Power" mode, similarly to what has 
been done successfully in the past, e.g. with ISOPHOT 
at 170 jtm on some sel ected low surface brightness fields 
llLagache fc DoleLl200ll) . However, MIPS is by design well 
calibrated for small spatial scales, e.g., point sources, be- 
cause the frequent stimulator flashes and the scanning 
strategy (acting like a chopping mode) properly track the 



fast r esponse of the photo conductors l|Cordon et all [2005. 
l2006t for instance). For these reasons, we do not use MIPS 
as an absolute photometer to estimate the level of the CIB 
and the foregrounds, but as a detector of small scale fluc- 
tuations to resolve the CIB based on 24 /im source obser- 
vations. Our approach has the important advantage that 
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it is not biased by the foregrounds and their modeling, 
which can lead to significant errors. 

We therefore start by reviewing the direct measure- 
ments, using absolute photometry in large beams, pro- 
vided mainly by the COBE FIRAS and DIRBE experi- 
ments and also the IRTS and rocket experiments in the 
near infrared (< 3/im). These measurements can be com- 
bined with indirect upper limits derived from observations 
of gamma rays from distant Blazars at TeV energies. 

To use the FIRAS and DIRBE data to provide CIB 
absolute measurements requires an accurate component 
separation. Local extended emission from interplanetary 
and interstellar dust can be removed using their specific 
SEDs and anisotropic spatial distributions traced indepen- 
dently, as well as time variability for th e zodiacal emission 
and scattering l|Hauser fc Dwekl l200ll for instance). Early 
gamma ray data from Blazars from the CAT experiment 
led to upper limits on the CIB intensity significantly lower 
than t he D IRBE re s iduals as pointed out bv lRenault et alJ 
l|200l|) a ndlWrightJ ll2004l). Recent results on more d istant 
Blazars l|Schroedteil 120051 lAharonian et all l2005|) con- 
strain the CIB even more in the near and thermal in- 
frared. Together with lower limits obtained by integrating 
the galaxy counts from HST, ISO, and Spitzer, these mea- 
surements tightly constrain the Extragalactic Background 
Light between ~ 0.8 to ~ 20 /im. 

At 160 /im, the CIB can be interpol ated from the 
DIRBE/COBE measurements at 1 00 /im j Lagach e et all 
|2000|) and 140 and 240 //m l|Hauser et all ll99S|) : 
0.78±0.21, 1.17±0.32, 1.09±0.20 MJy/sr, respectively. If 
the FIRAS photometric scale is used in the calibration 
(rather than the DIRBE photometric calibration), lower 
values are obt ained at 140 and 24 /mi: 0.7 MJy/sr and 
1.02 MJy/sr l|Hauser et all ll99Sh . A large uncertainty 
in the determinations at 100 and 140 fjxa comes from 
the zodiacal emission removal, as is also true at 60 /im. 
The DIRBE zodia cal emission model was obtained by 
iKelsall et alJ l|l998|) relying on the variability with viewing 
geometry. Its accuracy can be estimated a posteriori using 
the residuals observed at wavelengths where the zodiacal 
emission is at a maxi mum (12 and 2 5 /mi ). The residual 
emission, obtained by lHauser et alJ l)l998|) . has in fact a 
spectrum very similar to the zodiacal one. The residuals 
are about 4.7 x 10~ 7 Wm _2 sr _1 at 12 and 25 fim, far 
above the up per limit derived by high -energy experiments 
likeH.E.S.S. (|Aharonian et all T2005). but not very much 
larger than the uncertainties of the IKelsall et all (fl998) 
zodiacal emission model. A conservative estimate of the 
amount of zodiacal emission not removed in this model 
at 12 and 25 /im is ther e fore a bout 4 x 10~ 7 Wm^sr 1 . 
Using the Kclsall et al. ( 1998) smooth high latitude zo- 
diacal cloud colors, the amount not removed at 100, 140 
and 240 fim translates to 0.30, 0.14, 0.045 MJy/sr, respec- 
tively. This reduces the CIB from 0.78 to 0.48 MJy/sr at 
100 (im, from 1.17 to 1.03 at 140 /im and from 1.09 to 1.05 
at 240 /mi. Adopting the FIRAS photometric scale gives at 
140 and 240 /im, 0.56 and 0.98 MJy/sr respectively. From 
the above discussion, we see that the CIB at 140 /im - 



Table 1. Contribution to the CIB of S24 > 60 /iJy galax- 
ies. Sources of uncertainty come from photometry, calibra- 
tion, and large scale structure. The CIB column gives the 
best estimate from the discussion in Sect. 14.11 



A 






CIB 


% CIB 


/im 


nW m~ sr -1 


MJy/sr 


MJy/sr 


resolved 


24 


2.16 ±0.34 


0.017 ± 15% 


0.022 


79 


70 


5.93 ±1.02 


0.138 ± 17% 


0.15 


92 


160 


10.70 ± 2.28 


0.571 ± 21% 


0.82 


69 



the closest in wavelength to the 160 /im MIPS bandpass - 
is still uncertain by a factor of about 2 because of the un- 
certainty in the zodiacal level. The DIRBE/FIRAS mea- 
surement of the CIB at 240 /im suffers less from zodiacal 
residuals and photometric calibration uncertainty. 

A firm uppe r limit of 0.3 MJy/sr a t 60 /im has 
been derived by iDwek fc KrennricrJ ((2005) using obser- 
vations of TeV gamma r ay em ission from distant AGNs. 
Mivillc-Dcschcncs ct al ] i2002^l uses a fluctuation analy- 
sis of IRAS maps to set an upper limit of 0.27 MJy/sr 
and give an estimate of 0.18 MJy/sr, on the assumption 
that the level of fluctuations-to-total intensity ratio is not 
strongly wavelength dependent. 

At 24 /im we use for the contribution of I R gal ax- 
ies to the CIB the estimate of IPapovich et al.1 l|2004|) of 
2.7jlo'7 nWm _2 sr _1 . This value comes from 1) integra- 
tion of the source counts down to 60 /iJy giving 1.9 ± 0.6 
sr 



nWm- 2 ™- 1 



2) extrapolation of the source counts to 

2 „_-i. 



lower fluxes, giving a cont ribution o f 0.8±^i nW m " sr 
and 3) upp er limits fromlStecker fc De Jagerl 1(199 7*) and 
from CAT fo.ena.nlt, et all 12001 h . 

The most constraining measurements and lower and 
upper limits on the Cosmic Optical Background (COB) 
and the CIB from 0.1 /im to 1 m m are all reported in 
Figure El The lLagache et all l(2004h model predicts a CIB 
at 240 /im of 0.98 MJy/sr, which is in very good agreement 
with the estimate from combined measurements discussed 
above. Furthermore this model agrees with the observa- 
tional constraints (e.g. number counts, CIB intensity and 
fluctuations). We can thus take the CIB values from this 
model as a reasonable interpolation between the better 
constrained CIB values at shorter and longer wavelengths: 
0.82 MJy/sr at 160 /im, and 0.15 MJy/sr at 70 /im. 

4.2. Contributions from MIR Sources with 

S24 > 60/uJy 

To estimate the contribution of MIR sources to the back- 
ground, we add up the brightnesses of all the S24 bins to 
get the integrated light at 24, 70 and 160 of all the resolved 
24 /im sources. Each S24 bin is corrected for incomplete- 
ness. The results are presented Table Q and Figure 
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Fig. 10. Cumulative brightness (left axis, in MJy/sr) and 
fraction of the background resolved (right axis) at 160, 
70 and 24 /jm (from top to bottom), as a function of the 
24 (im flux in the three combined MIPS GTO Fields. A 
completeness correction has been applied based on the 
incompleteness level of the 24 /im sample. The horizontal 
dashed line represents 50% of the CIB at each wavelength. 



shows the cumulative integrated light from galaxies in the 
FIR as a function of 6*24. For a sanity check, we obtain 
that at 24 microns the percentage of the CIB that is re- 
solved is 79%, which is in agreement with lPapovich et alJ 
l)2004|) within the error bars. At 70 and 160 micron we 
resolve 92% and 69% of the background, respectively. 

Half of the 24 /mi CIB is resolved by sources with 1S24 > 
190/iJy. In the FIR, half of the 70 /im CIB is resolved by 
24 /im sources brighter than S24 ~ 220/iJy, and half of 
the 160 /im CIB is resolved by 24 /im sources brighter 
than 5*24 ~ 130/iJy. This difference between 70 and 160 
suggests that the CIB at 160 /im is dominated by galaxies 
at slightly higher redshift than at 70 /im, a consequence of 
the spectral shape of LIRGs and ULIRGs or, equivalently, 
the effect of k-correction. This point is illustrat ed in figures 
5 and 6 of the review bv lLaeache et all l|2005l) . 

To put in perspective the problem of resolving the CIB 
and what the stacking analysis accomplishes, we plot in 
figure ITTI the new observed constraints on the extragalac- 
tic background SED. The fraction of the CIB resolved at 
MIPS wavelengths by unbiased surveys w as 79%, 20% and 
7% at respective l y 24, 70 and 160 /im 1 Paoovich et all 
l2004 iDole et all l2004a| . When using the present stack- 
ing analysis, this fraction rises to 92% and 69% at 70 and 
160 /im respectively, and is represented by the red lower 
limits (see also Table 0. 



Based purely on observations without modeling of 
galaxy SEDs, we find that most of the FIR background 
is resolved into MI R galaxies. This co nfirms the model- 
dependent result of lElbaz et al.l (|2002). This analysis is 
the first direct resolution of the CIB simultaneously in 
the MIR and the FIR. 

Moreover, we can now securely establish the physical 
parameters of the typical galaxies responsible for most of 
the CIB near its peak. Previous studies bas ed on ISO al- 
ready characterized the 15 (an population llFlores et all 
IT991 lElbaz fc CesarsU Eool iFranneanViini et all 1200.4 
for instance); see iLaeache et all ()2005|) for a review. 
Our 24 /tm sample is almost complete in flux (80% 
completeness down to S24 = 80/iJy and 50% at 
60/iJy), and the physical properties o f S?a > 80/iJy 



galaxies have been extensively studied l|Le Floc'h eij al, , 
2004 IPerez-Gonzalez et all 120051 ILe Floc'h et all I200E : 
Caouti et all 12006^ . These works, mainly targeting the 
CDFS field, show that 25-30% of the 24 /im galaxies lie 
at redshifts z > 1.5, and that the redshift distribution 
peaks around z ~ 1 (between 0.7 and 1.1). Assuming the 
CDFS is a representative field, the MIR and FIR CIB is 
thus mainly composed of galaxies with typical redshifts of 
unity, with a contribution from z > 1.5 galaxies. At these 
redshifts, the galaxies are mostly LIRGs with typical bolo- 
metric luminosities of about 3 x 10 11 L© (between 10 11 
and 10 12 £ Q ) forming about 50 M^yr' 1 (20-130). They 
have intermediate stellar masses of about 10 10 to 1 11 Mq 
[jPerez-Gonzalez et all 120051: ICaputi et all l2006|) . From 
this latter work we can also estimate the specific star for- 
mation rates of these galaxies to be between 0.1 and 1 
Gyr^ 1 . 



4.3. Mean colors of the galaxies contributing to the 
CIB 



Looki ng at the 24 /tm number counts of IPapovich et all 
( 2004) , one can see that the bulk of 24 /im CIB is mainly 
due to sources with 130 < 6*24 < 400 /tJy. We select three 
cuts in S24 to investigate the colors of the contributions to 
the CIB by different galaxy populations. In the following, 
redshifts come from ICaputi et all l|2006j) (see their figure 
5), and the relative contributions come from the integra- 
tion of the IPapovich et all l|2004h source counts and the 
lLagache et alJ l)2004^ model. The cuts are: 

• Above 400 /iJy: bright galaxies contributing about 25% 
to the 24 /tm CIB. The redshift distribution has a mean 
of 0.53 and a median of 0.44. 

• 130 < S24 < 400 /iJy: galaxies contributing the most to 
the 24 /im CIB, about 30%. The redshift distribution has 
a mean of 1.18 and a median of 1.03. 

• 60 < S24 < 130 /iJy: fainter galaxies with relatively low 
contributions to the 24 /*m CIB (about 15%). The redshift 
distribution has a mean of 1.27 and a median of 1.11. The 
sample becomes incomplete at 60 < S24 < 80/iJy, so the 
mean redshift may be underestimated. 
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Table 2. Mean observed colors in I v of S24 > 60 /xJy 
galaxies contributing to the CIB. 



S24 in /iJy 


160/70 


160 / 24 


70 / 24 


S24 > 400 


3.2 ±0.4 


29.7 ±3.8 


9.6 ±0.8 


130 < S 2 i < 400 


4.4 ±0.5 


40.6 ± 4.2 


9.4 ±0.9 


60 < 5 2 4 < 130 


5.3 ± 1.6 


32.7 ±6.8 


6.3 ± 1.1 


CIB a 


5.5 


38.0 


6.9 



a Data and Model; See Sect.ETTl 



The observed colors change systematically with S24, 
as can be seen in the bottom of figure [7| <SVo vs S24 shows 
a slope larger than one, when <Si6o vs S24 shows an av- 
erage slope of the order of unity. Table |3 gives the colors 
and their associated la uncertainties in the three bins. 
The 160/24 color is compatible with a constant of 33. The 
70/24 color increases with the flux from 6.3 to 9.6. Finally, 
the 160/70 color steadily decreases with flux. 

These colors can be interpreted as the SED of a LIRG 
being redshifted, since fainter 24 /xm sources lie at larger 
redshifts: the 160/70 ratio increases (with decreasing flux) 
because the peak of the big grains' FIR spectrum is shifted 
longwards of 160 /im. The color ratios involving the 24 /im 
band are less obvious to inter pret, since the Polyc yclic 
Aromatic Hydrocarbon (PAH) ijPueet fc Leger[[l989|) fea- 
tures (especially between 6.2 and 8.6 /im) and the silicate 
absorption feature are redshifted into and then out of this 
band. The 70/24 ratio evolution might have for its origin a 
mix of PAH (increasing the 24) and very small grains con- 
tinuum (decreasing the 70) being redshifted, that cancel 
out each other. 

If one wants to extrapolate the contribution of fainter 
(5*24 < 60 /iJy) MIR galaxies to the FIR CIB, a conserva- 
tive approach is to use a constant 160/24 and 70/24 color 
ratio for the unresolved population. To set these ratios, we 
take the colors from the faintest population (60 < S24 < 
130 /iJy); this faint population presumably has the clos- 
est characteristics to the unresolved one. We will there- 
fore use Ii6o/^24 = 32.7 ± 6.8, and /70//24 = 6.3 ± 1.1 
(from Table[21). Since the contribution to the CIB of these 
faint galaxies is modest (30% at most), the large uncer- 
tainties in these color ratios will not dominate the total 
background estimate. 

5. New Estimates of the Cosmic Far-Infrared 
Background 

5.1. New Lower Limits at 70 and 160 fim 

The present stacking analysis performed on detected 
galaxies S24 > 60 /xJy gives strong measured lower lim- 
its to the CIB due to galaxies at 70 and 160 fim, without 
requiring any modeling. To determine upper limits to the 



FIR CIB requires a different approach. There are many 
difficulties at 70 /im in extracting an accurate value of the 
CIB, mostly due t o the problems in the removal of the zo- 
diaca l component ijFinkbeiner et all l200Ct iRenault et all 
l200ll for instance). At 160 fim the CIB estimate is more 
robust, but still with a significant uncertainty (factor of 
~3, see Sect. Ell . 

Another way to get a good estimate of the FIR galaxy 
CIB brightness is to estimate the unresolved 24 fim back- 
ground fraction, use the 160/24 and 70/24 colors measured 
for the weakest sources, and then apply these colors to the 
unresolved part to get the 70 and 160 /im background es- 
timates. Thus, we extrapolate the colors of galaxies with 
524 < 60 /xJy using the colors of the 60 < S24 < 130 fiJy 
galaxies derived in the previo us section. To estimat e the 
unresolved 24 /im background. lPapovich et all 1 20041) used 
a simple extrapolation of the differential number counts. 
Since the slope of the counts below 100 fiJy is strongly de- 
creasing (—1.5 ±0.1 in dN/dS), the integral is dominated 
by the largest fluxes £24- The estimate is robust, unless a 
hypothetical faint population exists. The remaining unre- 
solved 24 /im background created by S24 < 60, fiJy sources 

(to be compared to 2.16 



is therefore 0.54 nWm 
nWm 



2 sr 1 for S24 > 60,/iJy sources). 



We derive the extrapolated FIR CIB level due to IR 
galaxies using: 



I\ 24 

vl v {\) = i/I„(24) x -4- x — 

-<24 A 



(2) 



The results of the extrapolation are presented in 
TableEl We obtain 7.1±1.0 and 13.4±1.7 nWrn^sr' 1 , at 
70 and 160 /xm respectively. Our new estimate, based on 
the integration of all the 24 /im IR galaxies, is in principle 
a lower limit because it does not account for any diffuse 
emission unrelated to the IR galaxies, nor for a small frac- 
tion of IR galaxies that might have been missed. Indeed, 
the extrapolation in color of the unresolved 24 /im popula- 
tion accounts for the faint-end of the luminosity function, 
but not for the hypothetical very high-redshift sources, 
or faint local galaxies with high FIR output, like a hypo- 
thetical population of elliptical galaxies with large 160/24 
colors. However, if this population exists, its contribution 
to the FIR background is constrained by the upper limits 
to be less than ^20%. 

Our estimate a t 70 fim is higher than the 
lLagache et alJ l)2004l) model estimate by 11%, and lower 
by about 13% at 160 /im. About 25% of the CIB bright- 
ness at 70 and 160 fim comes from faint MIR sources 
(S24 < 60 fiJy). Assuming our new FIR CIB values repre- 
sent the actual CIB values, we estimate that our stacking 
analysis of S 2 4 > 60 fiJy galaxies finally resolves 75-80% 
of the background at 70 and 160 fim. We also show that 
the population dominating the CIB is made of galaxies 
seen at 24 /im and their simplest extrapolation to lower 
fluxes. 
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Fig. 11. Extragalactic Background Light Spectral Energy Distribution from 0.1 /im to 1 mm, with new constraints 
from MIPS. Red arrows (lower limits) represent the fraction of the CIB resolved at 70 and 160 fxm using the stacking 
analysis for sources with > 60 /iJy. The blue square represents the contribution of all the unresolved 24 /jm sources 
(extrapolation from number counts), and the blue arrows represent the contribution of all 24 /im sources to the FIR 
background, using a simple color extrapolation for 70 and 160 /im (Sect. 0. See Figure for the other symbols. 



Table 3. Contributions of the 24 /im galaxies to the FIR 
CIB in nW m~ 2 sr _1 . For the S24 < 60 /zJy galaxies, a 
simple color extrapolation has been used, as described in 
Sect. 





24 /im 


70 /im 


160 /mi 


> 60 /J, Jy 


2.16 ± 0.26 


5.9 ±0.9 


10.7 ± 1.6 


< 60/iJy a 


0.54 


1.2 ±0.2 


2.6 ±0.5 


total CIB 6 


2.7 C 


7.1 ± 1.0 


13.4 ± 1.7 


CIB prior c 


2.7 C 


6.4 C 


15.4 C 



a Estimate using an extrapolation from 60 to /zJy. 

b CIB estimate due to IR galaxies. 

c Data and Model; See discussion Sect. 14. ll 



5.2. Spectral Energy Distribution of the Extragalactic 
Background 

In the near and mid-IR, upper and lower limits tightly con- 
strain the EBL SED: 1) with HST+Spitzer and H.E.S.S 
between 0.8 and 4 /xm, and 2) with ISO, Spitzer and CAT 
between 5 and 24 /im. In this range, the EBL SED is 
constrained to better than 50% (and to the 20% level in 
several wavelength ranges) . The EBL is now also well con- 
strained in the FIR; direct measurements of the diffuse 
emission and our new lower limits constrain the CIB SED 
to the 50% level. 

The permitted zone for the EBL SED is presented in 
Figure lT2l This zone is defined as the area between current 
upper and lower limits. In this zone, the COB brightness 
ranges from 19.5 to 35.5 nW m~ 2 sr -1 , and the CIB from 
24 to 27.5 nW m~ 2 sr _1 . The ratio COB/CIB thus ranges 
from 0.7 to 1.5. 

From these constraints, we may derive a conservative 
estimate of the EBL SED, that typically lies between 
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Fig. 12. Cosmic Optical Background and Cosmic Infrared 
Background due to galaxies permitted zone estimate 
(shaded area), using upper and lower values. See Figure^! 
for the other symbols. 
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Fig. 13. Our best Cosmic Optical Background (blue- 
shaded) and Cosmic Infrared Background (red-shaded) 
estimates. The gray-shaded area represents the region of 
overlap. See FigureElfor the other symbols. 



the upper and lower limits and that makes use of well 
kn own physical p r ocesse s. The CIB estimate, based on 
the lLagache et all l|2004l) model, agrees with the data and 
is strongly constrained in the MIR and the 240-400 /im 
range. It strongly decreases with increasing frequency be- 
low 8 /im because of the main PAH features at 6.2 to 
8.6 /im being redshifted. The COB estimate also decreases 
with increasing wavelength above 2 /im because of the old 
stellar population SED. This simple SEP b e havior is in 
agreement with the model of IPrimack et alJ lll999l) . Our 
reasonable guess is that the COB and CIB have equal 
contributions around 8 /im. 

FigurclTSlsliows our smooth EBL SED estimate (thick 
line), as well as our best estimate of the COB (blue 
shaded) and the CIB (red shaded). The overlap region 
where both COB and CIB contribute significantly and 
the resulting total EBL is shown as the gray-shaded area 
around 8 /im. We find that the brightness of the COB is 
23 nW m~ 2 sr" 1 , and 24 nW m~ 2 sr" 1 for the CIB. The 
ratio between the COB and CIB is thus of the order of 
unity for this EBL SED. 

Our results are in contradiction with IWrightJ {2004) 
who finds a COB/CIB ratio of 1.7, and values at least 50% 
higher than ours: 59 nW m~ 2 sr" 1 (CO B) and 34 nW m~ 2 
sr -1 (CIB). However, the lWrightl <)2004|) estimate came be- 
fore t he strong upper limits of H.E.S.S ijAharonian et all 
2005) below 4 /im. This limit puts the COB much closer to 
the integrated light from galaxy counts than to the diffuse 
measurements. From the galaxy counts and stacking anal- 
ysis (lower limits), and high-energy experiments (upper 
limits), the EBL is now very well constrained. In partic- 
ular, we can now securely state that the contributions to 
the EBL of faint diffuse emissions outside identified galaxy 
populations - too weak to be detected in current surveys, 
like population III stars relic emission, galaxy clusters, 
hypothetical faint IR galaxy populations - can represent 



only a small fraction of the integrated energy output in 
the universe. 

5.3. The Extragalactic Background vs the Cosmic 
Microwave Background 

It is interesting to update the contributions of the most 
intensive electromagnetic backgr ounds i n the uni verse, as 
has b een done for instance by IScottl l|2000|) or IWrightJ 
(200^, and we schematically represent these in Figureim 
Obviously, the Cosmic Microwave Background (CMB) 
dominates the universe's SED, and accounts for about 
960 nW m -2 sr" 1 . We showed that the CIB and COB 
each account for 23 and 24 nW m~ 2 sr -1 , respectively. 
With a total of 47 nW m -2 sr -1 in the optical and the 
Far-Infrared, the EBL represents about 5% of the bright- 
ness of the CMB. Taking into account the complete SED 
of the EBL will not change this picture, since the con- 
tributi ons to the total EBL brightness of the radio, UV, 
X-ray IIMushotzk v et al.(. [2000: Has inger et all 1200 lj) and 
7 ray l|Strong et all l2004|) extragalactic backgrounds are 
smaller b y one to three orders of magnitude than the COB 
and CIB (IScottl 1200(1, . 

The galaxy formation and evolution processes provide 
5% in brightness of the electromagnetic content of the 
Universe. Half of the energy comes in the form of starlight 
(COB) and half as dust-reprocessed starlight (CIB). The 
maximum of the power distribution is at ~ 1.3 /im for 
the COB and - 150 /im for the CIB (Fig. HU). There 
are therefore on average 115 infrared photons for 1 visible 
photon emitted in these processes. 

6. Conclusions 

Our key points and results for the resolution and charac- 
terization of the FIR CIB and the EBL are: 
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• A stacking analysis in three fields covering 0.85 square 
degrees including a sample of 19181 MIPS 24 /um sources 
with S24 > 60 /iJy lets us probe faint 70 and 160 /im galax- 
ies one order of magnitude below the confusion level and 
with a high signal-to-noise ratio. We take into account in 
our noise budget uncertainties coming from: photometry, 
calibration systematics, and large-scale structure. 

• 24 /im galaxies down to S24 = 60/iJy contribute 79%, 
92%, 69% of the CIB at respectively 24, 70 and 160 
(using 2.7, 6.4 and 15.4 nW m~ 2 sr^ 1 as the total CIB 
values at 24, 70 and 160 /im, respectively). This is the first 
direct measurement of the contribution of MIR-selected 
galaxies to the FIR background. 

• We derive the contributions to the CIB by flux density 
bin, and show good agreement between our stacking anal- 
ysis and the published source counts. This is a strong con- 
straint for models. Moreover, we show that the CIB will be 
mainly resolved at flux densities of about <SVo ~ 0.9 mJy 
and Sieo ~ 3 mJy at 70 and 160 /im, respectively. 

• We directly measure that the total CIB, peaking near 
150 u m, is largely resolved into MIR galaxies. Other 
works lIPerez-Gonzalez et al.U2005ULe Floc'h et~alll2005t 
ICaputi et all I2006L especially) show that these MIPS 
24 /im sources are ~ 3 x 10 11 L Q LIRGs distributed at 
redshifts z ~ 1, with stellar masses of about 3 x 10 10 to 
3 x 10 11 Mq and specific star formation rates in the range 
0.1 to 1 Gyr" 1 . 

• Using constant color ratios 160/24 and 70/24 for MIR 
galaxies fainter than 60 /iJy, we derive new conservative 
lower limits to the CIB at 70 and 160 /im including the 
faint IR galaxies undetected at 24 /im: 7.1 ±1.0 and 13. 4± 
1.7 nW m~ 2 sr _1 , re spectively. These new estimates agree 
within 13% with the lLagache et all 1 2004h model. 

• Using these new estimates for the 70 and 160/im CIB, 
we show that our stacking analysis down to S24 > 60 /iJy 
resolves >75% of the 70 and 160 /im CIB. 

• Upper limits from high-energy experiments and direct 
detections together with lower limits from galaxy counts 
and stacking analysis give strong constraints on the EBL 
SED. 

• We estimate the Extragalactic Background Light 
(EBL) Spectral Energy Distribution (SED) permitted 
zone (between lower and upper limits), and measure 
the optical background (COB) to be in the range 19.5- 
35.5 nW m~ 2 sr _1 , and the IR background (CIB) in the 
range 24 to 27.5 nW m~ 2 sr -1 . The ratio COB/CIB thus 
lies between 0.7 and 1.5. 

• We integrate our best estimate of the COB and the CIB, 
and obtain respectively 23 and 24 nW m~ 2 sr _1 ; We find 
a COB/CIB ratio close to unity. 

• The galaxy formation and evolution processes have pro- 
duced photons equivalent in brightness to 5% of the CMB, 
with equal amounts from direct starlight (COB) and from 
dust-reprocessed starlight (CIB). We compute that the 
EBL produces on average 115 infrared photons per visible 
photon. 
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Fig. 14. Schematic Spectral Energy Distributions of the 
most important (by intensity) backgrounds in the uni- 
verse, and their approximate brightness in nW m~ 2 sr _1 
written in the boxes. From right to left: the Cosmic 
Microwave Background (CMB), the Cosmic Infrared 
Background (CIB) and the Cosmic Optical Background 
(COB). 
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